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I.  INTRODUCTION 
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Recently,  the  development  of  a new  automatic  gyrocoinpassing  tech- 
nique for  Pershing  II  (PII) , an  improved  version  of  the  Pershing  missile 
has  been  reported  [1-JfljJ.  The  goal  is  to  have  a faster  and  more 
accurate  azimuth  self -alignment  with  minimum  effort  from  human  operators. 

Among  the  reported  subjects  are  a new  gyrocompass ing  concept,  the 
theoretical  foundation,  implementation  scheme,  methods  of  on-line  data 
reduction,  computation  algorithms,  programming  considerations  for  soft- 
ware, laboratory  test  method  ane  result  evaluation,  and  techniques  for 
field  test  and  data  reduction. 

This  report  documents  the  result  of  several  error  analyses  for 
PII  inertial  measurement  unit  (IMO)  gyrocompassing  performed  during  the 
course  of  gyrocompassing  development.  The  analyses  have  been  beneficial 
for  the  development  as  follows: 

a)  The  process  of  analysis  revealed  more  insight  for  understanding 
the  underlying  principles  dictating  the  dependence  of  gyrocompassing 

error  on  various  source  errors.  1 

* i 

b)  The  analysis  yielded  a formula  for  estimating  gyrocompassing  j 

error  from  known  source  errors. 

i 

c)  The  analytic  result  obtained  provides  a background  for  devel- 
oping techniques  for  identification  of  source  errors  using  certain  test  j 

data.  1 


Gyrocompassing  is  performed  on  a g labelled  ZM7  before  missile 
launching.  The  operation  amounts  to  establishing  the  initial  alignment 
of  the  DflJ  platform  coordinates  with  respect  to  a geodetic  coordinate 
system.  There  are  many  error  sources  which  effect  gyrocompassing 
accuracy.  Among  them  are  sensor  anomalies,  ground  vibrations,  and 
measurement  noise.  Their  effects  will  all  be  analysed  bars.  To  make 
this  report  reasonably  self  contained,  the  gyrocompassing  s chasm  used 
will  first  be  reviewed. 

The  geodetic  reference  adopted  is  a triad  (V,S,A)  as  shown  in 
Figure  1 where  three  reference  axes  are  pointed  northward,  eastward, 
and  vertically  downward.  Ideally,  the  x,  y,  and  a axes  Of  the  IMP* 
platform  are  in  the  downrange,  right-hand  crossrange,  amd  d mansard 
directions,  with  x and  y axes  in  the  horizontal  plana.  Qymrnwp  easing 
is  performed  to  determine  the  azimuth  of  the  x axis  aieaaatad  with 
respect  to  the  geodetic  north. 

II.  THE  QYROCOMPA88ING  SCHBUE 

The  gyrocompassing  scheme  [1,5,7]  adopted  Is  referred  to  as  "two- 
position  offset  zero-torquing  gyrocompassing."  In  the  scheme,  platform 


fine  alignment  is  done  analytically  for  two  platform  positions  which  are 
90*  apart  in  the  horizontal  plane.  Alignment  in  the  first  position  is 
primarily  for  determining  the  drift  of  an  equivalent  north  gyro  which 
is  a fictitious  gyro  representing  the  combined  effects  of  x and  y gyros 
in  the  north  direction.  This  equivalent  north  gyro  becomes  an  equiv- 
alent east  gyro  in  the  second  position.  Alignment  data  obtained  in  the 
second  position  together  with  the  east  gyro  drift  allow  the  determination 
of  azimuth  heading  of  the  platform's  x-axis . The  word  "offset"  refers 
to  the  feature  which  permits  gyrocompasslng  with  the  platform  x-axis 
offset  from  north  at  any  azimuth  heading.  Zero  torqulng  refers  to  the 
concept  of  opening  the  leveling  loops  of  the  platform  during  data  taking, 
thus  reducing  the  effect  of  gyro  torquer  scale  factor  uncertainties.  A 
summary  of  the  gyrocompasslng  procedure  is  given  in  the  following 
paragraphs. 

Step  1:  In  the  first  position,  the  platform  is  approximately  level, 
with  its  x-axis  at  any  offset  angle  a with  respect  to  north.  A crude 
estimate  of  o is  made  using  any  beat  available  true  heading  (BATH) 
technique . 


I 


Step  2:  Both  of  the  platform  leveling  loops  are  opened  at  outputs 
of  two  accelerometers,  which  removes  the  torqulng  of  two  level  axes  of 
the  platform.  Torqulng  of  the  vertical  axis  of  the  platform  is  main- 
tained at  a rate  equal  to  the  vertical  component  of  earth  rate.  Outputs 
of  x and  y accelerometers,  which  are  in  the  fora  of  incremental  veloci- 
ties AV  and  AV  , are  measured, 
x y 

Step  3:  A rotational  transformation  on  AV^  and  AVy  is  applied  to 
give  and  AVg,  the  equivalent  incremental  velocities  in  the  north,  and 
east  directions,  respectively.  The  transformation  is  represented  by 


AV„  “ AV  cos  a - AV  sin  o 
N x y 

AV_  “ AV  sin  a + AV  cos  a 
Ex  y 


Cl) 

(2) 


Step  4:  By  integrating  accelerometer  outputs,  apparent  north 
velocity  V and  east  velocity  V , caused  by  the  tilt  of  platform,  a 

obtained.  They  are  related  to  platform  parameters  as  follows: 

'll  ■ W ♦ WJ  V2  - 1/6  “so  V3 

¥e  ■ -W  - l'2  “io'2-  1/6  deo  V3 


(3) 


(4) 


where  the  unit  for  velocity  is  in  g-seconds  and  g is  the  gravitational- 
acceleration. 

In  Equations  (3)  and  (4),  0^  and  8g0  are  the  initial  misalignments 

of  the  platform  about  the  north  and  east  directions,  respectively; 

and  Dg0  are  initial  platform  drifts  about  north  and  east  directions; 

ftA  is  the  vertical  component  of  earth  rate  ft,  and  t is  the  time  variable. 

Vg  and  Vg  are  obtained  at  1250  instants  during  a 240-sec  period, 

resulting  in  2500  data  pieces  from  which  platform  parameters  are  deter- 
mined using  one  of  two  available  least-square  regression  algorithms. 


Step  5:  A first  aslmuth  misalignment  t(>  is  computed  using 


MO 

♦ - - - — - - tan  L , 

T ft  cos  L HO 

and  the  north  gyro  drift  DR  using 

V * *io  - <l  - •“  ♦>  - 

In  Equation  (5),  L is  the  latitude  of  launch  site.. 
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(5) 
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Seep  6:  Both  leveling  loops  are  closed  and  the  I MU  platform  Is 
slewed  approximately  90°  about  its  vertical  axis  so  that  the  original 
platform-north  becomes  the  new  platform-east  and  the  original  platform- 
east  becomes  the  new  platform-south. 

Step  7:  Both  leveling  loops  are  opened  again  at  outputs  of  x and  y 
accelerometers.  A new  offset  angle  a is  estimated  by  BATH  for  this 
second  posltic.  Repeat  Steps  2,  3,  and  A for  a new  set  of  platform 
parameters. 

Step  8:  The  refi.ed  initial  azimuth  misalignment  angle  6^  is 
computed  using 


e - 

AO 


^PE0^2  ~ (1Vl 

0 cos  L 


- <0NO)2  tan  L 


(7) 


where  ( )^  and  ( >2  indicate  results  obtained  at  the  first  and  second 
positions , respectively. 

Step  9:  The  azimuth  misalignment  is  the  finely  determined 

angular  correction  for  the  crudely  estimated  offset  angle  a.  Therefore, 
the  initial  azimuth  heading,  H^,  of  the  platform's  x-axis  can  be  com- 
puted from 


H0  * » * 9 AO 


(8) 


The  final  result  of  gyrocompassing  consists  of  accurately  deter- 
mined initial  azimuth  heading  HQ,  initial  platform  misalignment  angles 

and  initial  drift  coefficients  and  D^.  These 


«HO  “■>  <W 


parameters  will  all  be  used  for  the  initialization  of  the  missile's 
navigation  and  guidance  system. 

III.  ERROR  SOURCES 

Various  error  sources  may  be  grouped  into  three  different  cate- 
gories. The  first  category  consists  of  all  error  sources  which  result 
in  heading-sensitive  gyrocompassing  errors.  These  error  sources  are  as 
follows: 


a) 

b) 


B and  the  variable  part  AB  , AB  , 
z r x y 


Accelerometer  scale  factor  uncertainties  K , K , and  K . 

* y ® 


The  fixed  part  of  accelerometer  bias  uncertainties  B , B , 

x y 

ABz  representing  the  differences  of 
tar  biases  at  the  first  and  second  gyrocompassing  positions. 


* V*  s':*  ^ 


c)  Gyro  heading-sensitive  drift  uncertainties  H , H , H . 

87  88 

■ d)  Y- accelerometer  non-orthogonality  4 which  is  in  the  xy- 

*y 


plane. 

a)  Gyro  input-axis  g-sensitive  drift  uncertainties  D^x  and  D^. 

f)  Gyro  spin-axis  sensitive  drift  uncertainties  and  D^. 

g)  Gravity  uncertainty  eg. 

h)  Latitude  uncertainty  e^. 

The  subscripts  x,  y,  and  z denote  the  respective  platform  axes  to  which 
sensors  are  associated.  Error  sources  in  this  category  are  assumed 
stationary  during  a short  time  period,  for  example , one  hour.  It  should 
be  noted  that  all  known  anomalies  are  not  considered  errors  because 
their  effect  can  be  compensated  by  software  during  data  processing. 

The  second  category  consists  of  all  sources  which  cause  base  notion. 
These  sources  Include  ground  vibrations  and  wind  gusts  which  are,  in 
general,  time- varying. 

The  third  category  of  error  sources  includes  measurement  noise 
which  may  be  considered  as  independent  white  noise.  Gyrocompasslng 
error  due  to  measurement  noise  has  been  analyzed  and  recently  reported 
by  the  authors  [1,3,7].  Therefore,  this  case  will  not  be  discussed 
further  in  this  report. 

IV.  HEADING  SENSITIVE  ERRORS 

A.  Accelerometer  Outputs 

Three  rectangular  coordinate  frames  will  be  used  as  follows: 

1)  The  geodetic  frame  with  axes  N,  E,  and  A as  shown  in 

Figure  1. 

2)  The  ideal  platform  frame  with  axes  X,  Y,  and  Z which  are 
pointed  in  the  downrange,  right-hand  crossrange,  and  vertically  down- 
ward directions. 

3)  The  actual  platform  frame  x,  y,  and  z which  is  misaligned 
from  X,  T,  and  Z respectively. 

The  misalignment  of  the  platform  is  represented  by  three  Euler  angles 
*x*  V*  11114  0*  rot*t#d  aboht  the  x,  y,  and  a axes,  respectively,  in 

that  order. 


The  transformation  between  (N,E,A)  and  (X,Y,Z)  for  an  offset  angle 
a,  as  shown  in  Figure  2,  is  given  by 


Figure  2.  Offset  a and  non-orthogonality  5 
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Y - T 
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where 


[cos  a sin  a 0 1 
-sin  a cos  a 0 I 

0 0 ij 


Transformation  from  the  ideal  platform  frame  (X,Y,Z)  to  the  actual  plat 
form  frame  (x,y,z)  represents  the  effect  of  misalignment  aaglea  6^,  0y, 

0 , and  the  non-orthogonality  angle  6 as  shown  in  Figure  2.  They  are 


I - T’  12  i 
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T3-  -S6xyC6xy°  ‘ <13> 

0 0 ij 

The  notations  C6  » cos  e and  S6  » sin  e are  used  to  save  space. 

During  prelaunch  gyrocompass ing,  the  missile  is  stationary,  and  in 
the  absence  of  base  vibratory  inputs,  all  acceleration  components  are 
due  to  the  earth's  gravitational  acceleration  g.  The  measured  accelera- 
tions along  (x,y,z)  frame  are  given  by 


T3T2Il|aK 


where 


vi  r° 

•e  - 0 


The  unit  for  acceleration  Is  the  number  of  g's. 
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By  modifying  Equations  (14}  and  (15),  effects  of  accelerometer 
scale  factor  uncertainties,  accelerometer  bias  uncertainties,  and 
gravity  uncertainty  can  be  taken  into  account.  Let  K be  a scale  factor 
matrix;  that  is 


K 0 0 

x 


0 0 K 


Considering  the  previously  mentioned  uncertainties.  Equations  (14)  and 
(15)  give 


a - [I  + K]  T 

y 


vhere 

T - T3  T2  Tx  . (18) 

The  bias  uncertainties  have  the  same  unit  as  acceleration,  which  is  the 
number  of  g's. 


Equation  (17),  with  its  exact  details,  is  rather  complex  and  hard 
to  analyze.  Fortunately,  in  practice,  angular  misalignments  are  suffi- 
ciently small  so  that  small  angle  approximations  are  valid,  namely, 
cos  0*1  and  sin  0=6.  Applying  the  approximation  to  angles  6 , 0 , 0 , 

* y z 

and  6^  in  Equation  (18)  and  neglecting  second  and  higher  order  terms 
result  in 
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Substituting  Equation  (19)  Into  Equation  (17)  and  expressing  the 
result  In  expanded  font  .yield 

fir  11"^ 

v :®*VS*’V!  V 

a I - (1  + K ) 6 (-1  -e  )j  +j  * ! 

; y ; y x e . y 

*,!  *« 

.11  iii 

The  third  order  variations  K 6 c and  K 8 e in  Equation  (20)  are 

y x 8 x y g 

neglected 


a-(l  + K+  e)8+B  ^ 

% * g y x I 

a - -(1  + K + e ) 8 + B l 

y y 8 * y r 

a ■ -(1  + K + e ) + B 1 

* * 8 * / 


By  examining  Equation  (21) , it  is  seen  that  6 , the  y-axis  non- 

orthogonality,  is  absent.  The  reason  is  that  the  effect  of  6 on 

accelerometer  outputs  is  second  order  in  nature.  Physically,  if  the 

platform  has  no  misalignment,  6 has  no  effect  on  accelerometer  out- 

xy 

puts  during  gyrocompasslng . 

It  is  desirable  to  retain  6 in  the  model  for  accelerometer  out- 

xy 

puts  so  that  its  effect  on  gyrocompasslng  can  be  monitored.  Instead  of 
dropping  all  second  and  higher  order  terms  in  obtaining  Equation  (19), 
those  second  order  terms  containing  6^  will  be  retained.  Therefore, 

Equation  (18)  gives 


Co  - 8t  So 
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The  desired  expression  for  accelerometer  outputs  is  obtained  by  substi- 
tuting Equation  (22)  into  Equation  (17) . 

s * » ♦ V V V '*  \ 


«,  - -a  ♦ ^ v ♦ », 

at  * -a  ♦ Rt  4 tg)  ♦ 


The  effect  of  6 appears  only  in  the  y-acceleration.  This  effect 
xy 

vanishes  when  the  platform  is  level.  The  outputs  of  accelerometers  are 
time-varying  because  misalignments  are  changing  due  to  platform  drift. 
The  dynamic  property  of  the  platform  will  be  analyzed  next. 

B.  Drift  Dynamics  of  the  Platform 

The  platform  of  a giaballed  IMU  is  tightly  slaved  to  its 
gyros;  therefore,  g\  o drifts  completely  appear  as  part  of  the  plat- 
form drift.  Let  D',  ' and  D'  denote  drifts  of  the  platform  about  its 

X z 

x,  y,  and  z axes,  respectively.  When  leveling  loops  of  the  platform 
are  open,  the  drift  of  the  platform  consists  of  not  only  gyro  drifts, 
but  also  earth  rotation  and  the  effect  of  platform  misalignment; 
therefore. 


fD’l 
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L zJ 

,Dz. 

A. 

where  R^  * R cos  L and  R^  ” -R  sin  L.  The  first  right-hand  term  of 

this  equation  represents  gyro  drifts;  the  second  term  represents  the 
combined  effect  of  earth  rotation  and  platform  misalignment.  The  trans- 
formation T2  in  Equation  (24)  is  given  by 


T 


2 


L y * j 

which  is  obtained  by  applying  small  angle  approximations  to 
Equation  (12). 

It  is  noted  that 


(25) 
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(26) 


•here  the  subscript  No"  implies  initial  value.  Using  Equation  (26)  in 
Equation  (25)  and  substituting  the  result  Into  Equation  (24)  give 
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are  initial  values  for  D’,  D*  , and  D', 

x’  y s 

Taking  the  Laplace  transform  of  Equation  (27)  and  expressing  the 
result  in  matrix  form,  result  in 
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The  characteristic  polynomlnal  of  Equation  (29)  la 

|.i  - - .3  * (02  . - #(.2  ♦ »2>  . 


The  solution  for  Equation  (29)  is 


ro’  <s)i 

ro*  ’ 

X 

xo 

D'  (a) 

- [si  - E]-1 

D' 

y 

t 

yo 

<8>. 

' 

D* 

L zo  J 

In  the  expanded  form. 
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+ 82  J 


‘] 


s 


Co 


CL  SL  So 


d;  (s> 


• D'  E 

xo 


CL  SL  Co 


2 

CL  SL 


, Co  + So"l 

a2+*2 J 


+ D' 
yo 


CL  SL  So  a - 
a2  + Q 


Co 


CL  SL  So 


+ D' 
zo 


[s2L  (1  - S2L)s] 

L-  a2  + n2J 
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By  taking  th*  inveraa  Laplace  t rant fore  of  Equations  (33) , (34) , and 
(35),  tine  response  of  drifts  is  obtained  as  follows: 

D* (t)  - [C2L  C2a  + (1  - C2L  C2a)  C(Gt))  D' 

X 

♦ (C^L  Ca  Sa  C(Gt)  ♦ SL  S(Ot)  - C?l  Co  Sal  D' 

70 

+ (CL  SL  Ca  - CL  SL  Ca  C<Qt)  - CL  Sa  S(Gt)J 
D* (t)  - (C2L  Ca  Sa  C (Gt)  - SL  C(Gt)  - C2L  Ca  Sa)  D’ 

y XD 

+ (C2!  S2a  ♦ (1  - C2L  S2a)  C(Gt) ) 

+ (CL  SL  Sa  C(flt)  + CL  Ca  S(Gt)  - CL  SL  Sa)  D*o 

D*(t)  - (CL  SL  Ca  - CL  SL  Ca  C(Gt)  - CL  Sa  S(ftt)l  D* 

z xo 

(CL  SL  Sa  C(Gt)  - CL  Ca  S(Gt)  - CL  SL  Sa] 

+ (s2l  + a - s2l>  c(ot)i  d;o  . 

The  data  taking  period  for  gyroconpasslng  is  sufficiently  short  so  that 
C(Gt)  * 1 and  S(Gt)  * Gt.  Therefore,  Equations  (36),  (37),  and  (38)  can 
be  adequately  approxiaated  by 

D'(t)  • D'  + D'  G.t  + D*  So  Q^t 
x xo  yo  A to  h 

D’(t)  - D'  - D'  fl,t  + D'  Ca  £L.t  » . (39) 

y ' yo  xo  A *o  H 

D’ (t)  - D'  - D’  Sa  G„t  - D*  Ca  G^t 
x xo  xo  h yo  n 

By  Integrating  Equation  (39)  the  tine  response  of  platfom  niaalignaent 
is  obtained: 


" 8*0  + “xo1  + 7 <v  “a  + Bw  ••  V 

y«  ■ V + V + T <■>»  «•  °» - D»  V 

v»>  • •„  + ®;0‘  - 4 *•  s ♦ ^ °*  v 


C.  Errors  in  Velocity 


Substituting  Equation  (40)  into  Equation  (23)  gives  accelerom- 

*y 


eter  outputs  which  contain  Kx,  Ky,  Kz,  B^,  By,  Bz,  eg , and  as 
follows : 


a » (1  + K + e 
x x g 


) fe  + D’  t + V Ca  n„  - D’  fl.M+B 
L y°  y°  2 ' zb  N xo  A J X 

[<6,« + v v + <°» + d;»  v 

+ 4 <Dy°  “a  + “io  S«  “g  + d;„  C“  r*N  S*y  ' Bio  »*  V] 


(41) 


a “ -(1  + K + e ) 

y y g 


+ B 


, - -(1  + K + e ) + B 
z * g 


(42) 


(43) 


(44) 


The  velocity  is  obtained  by  integrating  Equations  (41),  (42),  and  (43) 
as  follows: 

Vx  ' [B,  + * Kx  + £g>  9yo]  ' + [(1  ♦ + ‘g>  D;»]  4 

+ [w  + K,  + V c“  S - “io  V]  T 

vy'[*y-  <1  + W <8»,  + V V]* 

' [<*-  + Ky  + eg>  + Dyo  V]T 

' [(1  + Ky  + eg>  <Dyo  “a  + Ko  s»  °h  + D»o  c“.  S 4xy 

- V,v]t 

»,  - [»,-<»  + *,  + y)t  • 


(45) 

(46) 


Equations  (44),  (45),  and  (46)  give  velocities  along  three  plat- 
form axes.  The  velocities  derived  froa  a set  of  equivalent  platforn 
north  and  east  accelerometers  can  he  obtained  from  Equations  (44)  and 
(45)  via  a coordinate  transformation.  They  are 


H 


- V Ca  - V So 

* y 2 

* c“  ([*« ♦ <l  + *, + v v]* + [(1  * ** + V S.])  t 

t T/«  4 * \ /tvf  ' rv  _ Kf  A Clt3  1 


[“♦K,+6,)  <o;«  ^ v]t  ) 

- *•  [»T  - « * S + V <*»  + V V]  1 

- « + S + •*>  <D^  + Djo  V]  T 

- (1  + K e ) (O’  ft.  + D’  Sa  ft*.  + D’  Ca  ft„  6 

' y g ' yo  A so  M so  n xy 


- D‘ 
xo 


V*  - V 8a  4-  V Ca 

* * y 


+ Ca 


D 


B * (1  + K + eB)  (6  +0  6 )li 

y ' y 8 xo  yo  xy  J 


- Ta  + K + e ) <D’  + D*  6 )1  V 

L y 8 *o  yo  xy  J 2 

- Td  + K + e ) (D’  ft.  + D’  Sa  ft*.  + D’  Ca  ft_  6 

J;  y g yo  A so  71  so  H xy 

- D'  ft.  )1  %-i  . 

xo  A xy  j 6 J 


(47) 


j 

2 

s»  {[*« + (i  + K. + v v]  * + [(1  + ** + v ■>;<.]  t 

+ [«  + K, + V S ' D«  v]  r) 


(48) 


Hie  error-free  worth  end  esst  velocity  can  be  obtained  froa 

Equations  (47)  and  (48)  by  nulling  l , 1 , B , B , e , and  d . In 

* y * y • ” 

addition*  D'  is  set  to  sero,  since  this  Is  one  of  the  assumptions  andar 
so  j 

which  gyrocoapassing  Equation  (7)  la  derived.  Denoting  error-free 

velocities  by  superscript  "0",  the  following  is  obtained:. 


£7 


* *1  . A . 


v„  - Ca  \ 3 t + D*  -=r-  - D’ 


r D n.  i 

yo  2 xo  A 6 I 


+ s“(v  + D;T  + D;oaAT 

“ (eE0t  + °E0  T ~ dno  sat| 


I t2  t3 1 

Sa(eyot  + D;oV-D^^Vj 

I c2  t3 

Cal  -e  t - D’  V - D'  tlA  =7- 
ixo  xo  2 yo  A 6 

" °N0  T “ deo  nA  T J 


where 


0 ■ 0 Ca  ~ 0 Sa  / c i \ 

NO  xo  yo  A31A 

6E0  " exo  50  + 9yo  Ca  (52) 

°N0  " °xo  Ca  * Dio  So  (53) 

°E0  “ Dxo  Sa  + Dyo  Ca  * (54) 

Equations  (49)  and  (50)  are  in  fact  Equations  (3)  and  (4),  respectively. 

Expressions  for  velocity  error  and  AV_  are  obtained  by 

0 o * 6 

subtracting  from  and  Vg  from  Vg  using  Equations  (47)  through  (50) 
as  follows: 

AVn  - Ca  - By  SaJ  c + jcK^  - Ky)  Sa  Ca  + Qg  S2a  6 J t2 


+ Id  - 0 

I X xo 


s»  - •„  S to]  % T 


Vg  - |bx  Sa  + By  Cajt 

+ k s »2» + °h  c2» + c,  % - °* 


Sa  Ca  5 


li! 
*J  2 


The  lest  bracketed  ten  of  Equation  (55)  comes  from  D*0  using  the  rela- 
tionship of  Equation  (28c)  but  with  known  quantity  deleted*  Conquering 

Equations  (55)  and  (56)  to  Equations  (3)  and  (4) , it  can  be  seen  that 
errors  in  the  determination  of  8^,  8^,  D^,  and  D’^are  given  by 


Ae*0  * *<BX  Sa  8y  Ca)  “ -*E 
A6BD  " *x  Co  - Ey  Sa  - B„ 
6B^>  - 6^  ^ Sa  Co  - kx  fl,, 


S2.  - Ey  n,  C2.  - eg  0, 


ADJo  - <5^  ^ S a + (Kx  - Ky)  Sa  Ca  . (40) 

D.  Effect  of  Heading-Sensitive  Gyro  Drift 

The  heading-sensitive  drifts  for  x and  y gyros  are 
Hsy  and  Hs»*  respectively.  Drifts  for  a set  of  equivalent  north  and 
east  gyros  are 


“s«  • V*  - Vs.  W1) 

»SE  • M«S“  + HayCa  ‘ «2> 

In  the  first  platform  position  of  gyrocompsssing,  the  determined 
drift  (Dg)^  of  the  equivalent  north  gyro  also  contains  the  heading- 

sensitive  drift  (H  ),  and  the  two  cannot  be  separated.  This  (D„). 

becomes  (Dg) 2 at  the  second  platform  position.  Because  of  the  heading- 

sensitivity,  (H  ).  is  not  equal  to  (H  ).,  although  both  are  referred  to 

the  same  physical  axis  of  the  platform.  The  difference  of  (H  )_  and 

SE  — - 

(H^i  is  contained  in  the  determined  (0^)^  at  the  second  position  and 

is  not  readily  separable  from  the  latter.  Fortunately,  separation  of 
(H$n) i from  (Djj)^  and  of  (Hsy>2  from  (Dy)2  is  not  necessary.  It  is  the 

relative  value  rather  than  the  absolute  value  of  the  heading-sensitive 
drift  which  will  be  used  in  the  gyrocompsssing  error  analysis.  Therefore, 
Equation  (61)  can  be  ignored  hut  Equation  (62)  should  he  Included  in 
Equation  (60)  to  give  the  error  in  the  determined  east-axle  drift  of 
the  platform  as 

**»»  • % s ♦ <*k  * V V*  *•  * 

- ^ “ V S *®  * *m  80  * **  ' <*** 


E.  Effect  of  Gyro  Input-Axis  g-Sensitive  Drifts 


The  PII  DAI  employs  a tuned-rotor  type,  two-degree-of-freedom 
gyro,  for  sensing  attitude  rate  about  the  platform's  x and  y axes.  From 
the  physical  consideration  of  the  device,  it  is  recognized  that  the 
coefficients  of  two  input-axis  g-sensitive  drifts,  Dj  and  D , are 


equal.  Hence,  both  of  them  will  be  denoted  by  D. 


Ixy* 


During  gyrocom- 


pass ing,  the  equivalent  north  axis  of  the  platform  has  little  tilting 
motion,  while  the  equivalent  east  axis  is  tilting  appreciably  due  to  the 
north  component  of  eaith  rotation.  The  tilt  of  the  east  axis  causes 
appreciable  g-sensitive  drift  for  an  apparent  east  gyro,  represented  by 
the  coefficient 


IE 


D_  So  + D_  Co 
lx  Iy 


DIxy  *Sa  + Co> 


(64) 


The  unit  for  all  D coefficients  is  degrees  per  hour  per  meter  per 
1 2 

second-squared  I (deg/hr) /(m/sec  )]. 


The  tilt  angle  for  the  platform's  east  axis  has  a nominal  value  of 
ftgt  when  leveling  loops  are  open.  The  average  drift  over  a time  Interval 

T,  the  data  taking  period,  is 


AIE‘  " T o °IE  8 V dt  “ 2 °IE  8 °NT 


(65) 


This  drift  should  also  be  included  as  part  of  the  error  of  the  determined 
east  axis  drift  of  the  platform  at  the  second  platform  position.  This 
is  done  by  including  Equation  (65)  in  Equation  (63)  giving 


ADEo“  6xy  S S ° + (Kx  ' V 


+ H<j[Sa  + H<yCa+i  DJxy  (So  + Co)  g J^T 


(66) 


where  Equation  (64)  has  been  used  to  express  D^g  in  terms  of  DIxy* 


F. 


Effect  of  Gyro  Spin-Axis  g-Sensltive  Drifts 

The  effect  of  spin-axis  g-sensitive  drifts  about  the  platform's 


x and  y axes  can  be  analyzed  similarly.  For  the  PII  IMU,  the  spin-axis 
of  level  gyro  is  oriented  vertically  and  experiences  an  acceleration  of 
1 g.  The  effect  of  the  slight  tilt  of  the  vertical  axis  may  be  Ignored 
because  of  the  insensitivity  of  the  cosine  function  near  null.  The 
characteristics  of  the  spin-axis  g-sensitive  drifts  about  the  x and  y 

The  equivalent 


axes  are  represented  by  coefficients  Dgx  and  D#y 


coefficients  for  the  platform's  north  and  east  axes  are  given  by 
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i i 


! i 


H 


AV 


t<%)>2  ' OVl1  a CL  + M£g>2  ' <VlJ  0 * cl 

a*c\ 


- *eH0  t4Kl  L " eM0  «*A  ^ 

Using  Equations  (6) , (53) , (28*) , and  (28b) , the  variation 
be  obtained  as  follows ; 

<“m)  • w"io>i  *“»«•-  % s“ 

* Ca  {-<«  Sa  - Ca)  11SL  + 9 OCU  t, 
so  yo  i> 

- Sa  {-(6  Ca  + Sa)  flSL  - 6 OCL}  e. 

*0  XO  L 

- (asL  ♦ qcl)  tL 

Similarly,  Equations  (54),  (28a),  and  (28b)  give 

AD'  - AD’  Sa  + AD*  Ca 
KO  so  yo 


- Sa 

Sa  - Ca)  QSL  + 6 QCL)  c, 
yo  L 

(75) 

i 

t 

♦ Ca  {-(S^ 

Ca  ♦ Sa)  QSL  - 0 OCL}  t. 

XO  L 

! 

• H os 

v so 

ocl)  cL 

(76) 

\ 

1 

J 

| 

It  is  noted  that 

i 

i 

“no-0 

(77) 

* 

with  respect  to  as 

can  base  seen  fron  Equation  (51). 

Substituting 

J 

Equations.  (75),  (78),  and  (77)  into  Equation  (73),  and  neglecting  all 
second-order  tarns  glee  tbs  error  of  asiauth  heading  In  a neat  equation 


A®go  * -eL  *»  1 


H.  The  Eaaultant  Gyroconpaaaing  Error 

The  effects  of  all  heading-  sensitive  errors  In  maccsfeislag 
mm  he  eegfci*ed.  haferring  to  Elgnre  3*  1st  the  asiauth  offset  is 
the  fleet  p let  Corn  position  he  a - 98s  and  in  the  ascend-,  poe4nian.be  a. 
8igce  cos  Co  - 9A*>  - aln  a and  ala  (a  - »*)  * -fin*  «.  Iqesjfana  (4), 
<m»  and  (71)  glee 


> .r  -•  « * . 

• ..  . * *ey  , 


>«>»•  v. 
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Figure  3.  Two  platform  positions. 

<4Vi  ’ <1Dm>i  - 4*  s*  - <4Wx  °*  “ <“ ioh  - (46eo>i  “a 

■ -5xy  «,  S“  <=“  - *x  °»  °2“  - Ky  °N  s2“ 

- eg  £!„  + 8 <Dgx  Sa  + Dgy  Ca)  - (Bx  Sa  + By  Co)  0^  (79) 

and  Bquation  (72)  gives 

(AD’  ),  « « JL.  S2o  + (K  - K ) So  Co  + H So  + H Co 
iyj  l xy  n x y w bx  ay 

+ 1/2  DJxy  (So  + Co)  gl^T  + g (D^  So  + Co)  . (80) 

In  Section  111,  it  was  pointed  out  that  accelerometer  bias  uncertainties 
might  be  different  between  the  first  and  second  gyro compassing  positions 
Therefore,  if  Equation  (57)  gives  AO^  for  the  first  position,  then  for 

the  second  position 


<A6wn>9  - “(B  + AB  ) So  + (B„  + AB  ) Co 


(81) 


1 


.‘The  resultant  gyroconipasslrig  error  i$  the  sun  of  the  error  due  to 
er  arid  errors  due  to  dll  ot'^ter  sources,  from  Equations  (7)  rind  (78), 

this  Is  given  by 


A0 


2 " (AVl 


AO 


+ feL  - (A«j^)2J  tan  i 


(82) 


Using  ( A6^j)  2 from  Equation  (81)  arid  A(D^)  j rihd  AXfw^  from 
Equations  (79)  arid  (80)  and  rioting  that  tan  L - yield 


A0irt  " K (So  Ca  + C2d)  + K (S2a  - Sa  Co)  + AB-  So  tan  L 
ao  x y * 


+ AB  Co  tari  L + H -v 
y sx 


So/Oj, 


+ ft 


Sy 


Co/Ojj  + (S  o + So  Co) 


+ 1/2  D 


Ixy  (So  + Co)  g T + eg  + eL  tan  L 


(83) 


This  Equation  does  not  contain  B , 8 , D , and  D.  terms,  indicating 

T y ®y 


that  in  a first-order  analysis,  the  azimuth  heading  obtained  by  two- 
position  gyrocompassing  is  independent  of  fixed  accelerometer  bias  uncer- 
tainties arid  sfiri-axis  g-serisitive  drifts. 


Equation  (83)  is  the  final  rfcsult  sought.  Although  the  object  was 
to  analyze  all  heading-sensitive  errors,  the  result  dhows  that  errors 


due  to  e_  and  are  not  heading-sensitive  as  indicated  in  Equation  (83). 


The  set  of  units  for  various  quantities  in  Equation  (83)  can  be  traced 
out  through  the  derivation.  The  units  are  as  follows: 


A9aq  in  radians 


and  Ky  in  g’s/g 


AB  and  AB_  in  g's 

x y 


H and  H in  degrees  per  hour 

sx  sy 


B in  radians 

*y 


D in  (degrees  per  hour)  per  (meters  per  second  ) 


iri 


of  g's 


in  rad  laris 


% 


in 


f*r 


1* 


“-  --'--i --  '*■—  a&Wi 1 


■ L- .. 


T in  hours 


r 

I 


a i 


r 

* 

V 


g in  meters  per  second  . 

A set  of  convenient  units  which  is  more  often  used  in  practice  is 
given  as  follows: 

A0._  in  arc  minutes  (arc  min) 

AO 

K and  K in  micron's  per  g (iig/g)  or  parts  per  million  (ppm) 

* y 

AB  and  AB  in  micro-g's  (wg) 
x y 

H and  H in  degrees  per  hour  (deg/hr) 
sx  sy 

6 in  arc  seconds  (arc  sec) 
xy 

in  degrees  per  hour  per  nominal  g 
e in  micro-g's  (pg) 

in  arc  seconds  (arc  sec) 

^ in  degrees  per  hour  (deg/hr) 

T in  seconds.. 

The  g factor  in  the  D term  is  replaced  by  k with  a unit  in 

lxy  g 2 

number  of  nominal  g's  where  the  nominal  g is  9.8  m/sec  . 

Using  the  new  set  of  units.  Equation  (83)  becomes 

A0 . n - k.  [K  (So  Ca  + Co2)  + K (So2  - So  Ca)J 
AO  lx  y 

+ k,  [AB  So  + AB  Ca]  tan  L + k,  [H  So  + H Co]  /SL. 
l x y 2 sx  sy  N 

+ k,  6 [So  Co  + So2]  + k.  D_  [So  + Co]  k T + k.  e 
3 xy  a lxy  g 1-8 

+ kj  tan  L . (84) 

Values  of  k^,  i ■ 1 to  4,  are  obtained  as  follows: 

k . 57.3  (deg/rad)  * 60  (arc  aln/deg)  _ 3>438  x 1Q-3 
1 106  (parts) 

k2  - 57.3  (dag/rad)  * 60  (arc  min/deg)  - 3438 


5Z.-3.J 

57.3 


k4  - 1/2  x 


7*? 


60  (mta/deg)  x 3606  (sec/hr) 


0.07799 


kg  is  9.8  m/sec  , 

I.  A Sample  Calculation 

To  illustrate  the  nature  of  the  heading  sensitivity  of  indi- 
vidual error  components  and  their  .contribution  to  the  resultant  heading 
error,  a sample  case  la  computed  using  the  following  set  of  parameter 
valves. 

K a K • 100  ppm 
x y rr 

ABx  - ABy  - 100  yg 

H - H - 0.001  deg/hr 
sx  sy  ^ 

6 » 10  arc  sec 

xy 

DIxy  " O'0015  deg/hr/g 

e - -10  yg 


eL  - -10  arc  sec 

k - 0.999665 
8 

Q » 15  deg/hr 
L - 34.6425" 

T - 240  sec  . 


Figure  4 shows  the  variation  of  individual  heading  errors  as  s function 
of  nominal  heading  angle  a.  Figure  5 shows  the  variation  of  the 
resultant  heading  error.  All  individual  errors  have  sinusoidal  varia- 
tions except  those  due  to  and  e, . Errors  due  to  e and  e,  have 

8 nu  Eh 

const ant  values  Indicating  their  insensitivity  to  heading.  Errors  due 

to  K a K , and 6 have  non-tero  naans,  while  those  due  to  AB  . AS  , H , 

*7  *7  * 7 ax 

■ * eadJBj  hone  aero  naans. 
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Figure  4.  Individual  heading  errors 


By  examining  Figure  5 it  is  seen  that  gyrocompassing  error  is 
smaller  for  certain  headings  than  for  others.  This  phenomenon  suggests 
that,  for  each  IMU,  a nominal  heading  position  a can  be  identified  which 
gives  the  best  gyrocompassing  accuracy. 

V.  ERROR  DUE  TO  BASE  MOTION 

Base  motion  is  a random  process  which  is  lew  frequency  in  nature 
as  compared  to  the  instr  ment  noise;  therefore,  it  is  colored.  Its 
effect  on  the  accuracy  of  gyrocompassing  depends  on  the  details  of  the 
data  processing  algorithm  used.  Experiments  with  two  different 
algorithms  are  presently  under  way.  The  first  one  involves  the  use  of 
a conventional  least-square  polynomial  fit  formula  to  fit  tvo  second - 
order  polynomials  [12].  The  second  algorithm  employs  a parameter- 
correlated  least-square  regression  to  fit  a set  of  two  third-order  poly- 
nomials [2, 3, 5, 6, 7].  Here  the  error  analysis  is  performed  based  on  the 
first  algorithm. 

Equations  (3)  and  (4)  are  seen  to  be  third-order  polynomials  in  t. 
The  first  algorithm  uses  the  following  approximation  at  the  outset 


V - a + bt  + ct  . (85) 

The  determination  of  heading  is  dominated  by  the  first  term  of  gyro- 
compassing Equation  (7).  Therefore,  the  accuracy  of  gyrocompassing 
primarily  relies  on  the  determination  of  drift  which  is  proportional  to 

2 

the  coefficient  c of  the  t -term  in  Equation  (85).  The  present  analysis 
will  be  concentrated  on  the  c coefficient. 

The  conventional  least-square  regression  for  a second-order  poly- 
nomial is  to  minimize 


V (t)  - (a  + bt  + ct2)]2  dt 


31  31  31 

Letting  — ■ ■ 0 and  performing  several  integrations  of  the 

polynomial  give 

r r 2 3ir  1 

I,  T T I-  . 

1 T TT  a 
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2 T 3.  4 h 
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. J L 3 T 5 ; 


Solving  Equation  (84)  for  the  coefficient  C gives 


r . SO  r 180  . . 18G 

c7ii7i!7ii 

Squaring  equation  (89)  gives 


. /3£\2  T 2 . /180V2  T 2 /160V2  , 2 

\t3)  1 (l4)  2 \7)  3 

10800  _ . . 10800  , , 64800  _ , 

- ix  i2  + —r  i2  - -—9  - 12  13 


<8! 
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which  will  be  need  later  for  Obtaining  the  standard  deviation  of  C.  11 
base  notion  error  will  be  analysed  by  first  using  a velocity  aodel  and 
then  an  acceleration  nodal. 


A.  Analysis  via  Velocity  Model 

It  is  assuned  that  the  velocity  of  base  notion  has  the 
following  aodel: 

II 

V * ^n-1  An  8ln  (V  + ♦n)  * (9 


This  aaounts  to  representing  the  velocity  by  its  Fourier  coaponeats. 
Using  Equation  (91)  in  Equation  (88)  yields 


*2  - 


♦n>) 


(91 


•*»  <v  * v ■ •**♦»*  V **■  ^ * v) 


(«* 


- I*  |2  cos  *n  - 2 cos  + V 

2 2 1 
+ 2 w t sin  (w  t + <J>  ) - to  t cos  (u  t + 4)  J 
n n n n n n / 


(92c) 


Squaring  Equation  (92a)  and  dropping  cross-frequency  product  term*  yield 


(%)["*  (“»t  + *»)]| 

« 1 ( — | ( cos2  4 + cos2  (to  t + $ ) 

n-I  \V  I n n n 

N /A  \2  | 1 + cos  2 4 1 + cos  2 (w  t + 4 ) 

■UWI— 


2 2 

cos  4 + cos  (to  t + 4 ) - 2 cos  4 cos  (ut  + * ) 
n n n n n n 


+ 2 cos  A cos  (to  t + A ) 
n n n 


The  cross-frequency  product  terms  are  dropped  because  their  average 
values  are  zero.  Again,  all  terms  in  Equation  (93)  which  have  a zero 
time-average  are  discarded: 


(94a) 


To  be  conservative,  let 


N /A 
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Similar  considerations  produce 


(94  b) 


When  Equation  (94)  is  substituted  into  Equation  (91),  several  te 
cancel  each  other.  The  result  is  the  time-averqge  of  : 


Assuming  ergodlclty,  time-average  equals  enseafele  average  aad 

2 

Equation  (95)  gives  oc  , the  variance  of  C.  Since  the  north  aad 
channels  of  the  platform  are  similar,  the  following  la  obtained: 


2 2 
°S  “0CE 


Coefficients  and  Cg  are  related  to  drifts  by 


- -2  C„ 


“EO^S 


(94e) 


<94f> 


(94  g) 


which  can  be  obtained  by  compering  Equation  (85)  to  Equations  (3)  and 
(4).  The  standard  deviation  of  the  heading  error  can  be  obtained  from 
Equations  (7)  and  (97)  as 


. , ’**  r»#« 

■*>  . ■ ,r-  . -• 


where  the  second  term  of  Equation  (7)  has  been  neglected.  Using 
Equation  (96)  in  Equation  (98) 


(99) 
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or 


! 
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2/3 

/30\ 

0 cos  L 

(t3) 

(57.3  x 60)  arc  min 


(100) 


Equation  (100)  shows  that  the  error  due  to  base  motion  decreases 
with  frequency  u>n  and  time  T over  which  least-squares  regression  is  done 

Also,  the  base  motion  error  is  larger  at  higher  latitude. 


B.  Analysis  via  Acceleration  Model 

The  acceleration  model  for  the  base  motion  is  assumed  to  be 


a 


N 

y g (u  t + ♦ ) 
L . n n Tn 
n-1 


(101) 


Following  a similar  derivation  as  that  for  the  velocity  model  yields 


(i)  ■ n « 
o o 

where  co  ■ it  rad/ sec  and  n - 1 to  100,  Let 
o 

Gq  - 2 x 10"3  g for  all  n 


T • 240  sec 


(105) 
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VI.  CONCLUSION 


There  are  three  categories  of  gyrocompassing  errors.  Analysis 
I for  the  first  two  categories  has  been  treated  in  this  report;  the  third 

' category  has  been  reported  elsewhere  by  the  authors.  The  analysis  has 

f been  beneficial  in  three  ways.  It  offered  an  opportunity  for  a better 

I understanding  of  the  underlying  principle  governing  each  error  source. 

The  analytic  expressions  obtained  enable  an  estimation  of  gyrocompassing 
error  with  known  error  sources . It  also  provided  a model  which  can  be 
used  to  develop  identification  techniques  for  error  sources  based  on 
* test  data. 

s'  Errors  of  the  first  category  include  all  those  which  are  short-term 

] stationary.  A deterministic  type  of  analytic  expression , Equations  (83) 

or  (84),  for  heading  error  was  obtained  for  this  category.  Using  this 
expression,  identification  of  error  sources  can  be  devised.  By  compen- 
sating for  the  identified  error  sources,  gyrocompassing  accuracy  can  be 
Improved.  A scheme  for  accomplishing  this  will  appear  in  a later  paper. 

Errors  of  the  second  category  result  from  base  motions  which  are 
random  and  time-varying  disturbances.  Their  effect  on  gyrocompassing 
accuracy  has  been  analyzed  statistically.  Analytic  expressions  in 
Equations  (100)  and  (104)  give  the  standard  deviation  of  heading  error 
! in  terms  of  spectral  properties  of  the  base  motion.  Error  sources  of 

this  type  cannot  be  compensated  by  adjusting  hardware  parameters.  How- 
1 ever,  their  effect  can  be  minimized  by  using  a proper  data  processing 

i algorithm.  The  low-frequency  nature  of  these  types  of  disturbances 

implies  that  they  are  not  white;  therefore,  they  are  more  difficult  to 
filter  than  white  noise. 

I 

The  third  category  of  error  sources  includes  all  instrument  noise 
which  usually  can  be  approximated  as  white  noise.  A statistical  analysis 
for  this  category  has  been  reported  elsewhere  recently  and  is  not 
'}  detailed  in  this  report.  The  effect  of  this  type  of  noise  on  gyrocom- 

passing accuracy  can  be  reduced  by  proper  filtering. 

Finally,  it  should  be  pointed  out  that  a fourth  category  of  error 
sources  has  not  been  mentioned  in  this  report.  This  category  includes 
all  software  errors  such  as  round-off  error,  truncation  error,  and 
errors  induced  by  approximations  made  during  algorithm  development. 

This  last  category  of  errors  will  be  treated  in  the  future. 
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